Electron-conformational transformations in nanoscopic RyR channels govern both the 

heart's contraction and beating 



O 
(N 



Oh 
I 

o 



C/3 

o 

• I— I 

Oh' 



> 

O 

rn 
O 



X 



A.S. Moskvin,! A.M. Ryvkin,^ O.E. Solovyova,^.^ and V.S. Markhasin^ 

'Department of Theoretical Physics, Ural State University, 620083 Ekaterinburg, Russia 
^Institute of Immunology and Physiology, Ural Branch of RAS, 62004-9, Ekaterinburg, Russia 
^Department of Computational Mathematics, Ural State University, 620083 Ekaterinburg, Russia 

(Dated: January 18, 2013) 

We show that a simple biophysically based electron-conformational model of RyR channel is able 
to explain and describe on equal footing the oscillatory regime of the heart's cell release unit both in 
sinoatrial node (pacemaker) cells under normal physiological conditions and in ventricular myocytes 
under Ca^"*" SR overload. 



Calcium (Ca^+) dynamics is of a principal impor- 
tance for functioning of different heart's cells from atrial 
and ventricular cardiomyocites to sinoatrial node cells 
(SANC) though the former are responsible for the heart's 
contraction while the latter for primary heart's pacemak- 
ing, respectively—. Cardiac contraction in cardiomyocites 
is activated by an increase in intracellular calcium con- 
centration (Ca^"*"), most of which comes from a specific 
calcium cistern of sarcoplasmic reticulum (SR). Ca^"*" is 
released via the ryanodine receptors (RyR) in response 
to Ca^"*" entering the cell via the L-type channels (see 
Fig. 1). The cardiac type RyR is the common major Ca^"*" 
release channel type in SANC and ventricular myocytes. 
It has been experimentally documented in chemically 
skinned and voltage-clamped SANC, in which effects of 
volt age- activated sarcolemmal ion currents are excluded, 
that the isolated SR is capable to spontaneously and 
rhythmically release Ca^"*" via RyRs-Si^,. These spon- 
taneous, rhythmic, local subsarcolemmal Ca^"*" releases 
(Ca^"*" clock), which occur in SANCs, interact somehow 
with the classic sarcolemmal voltage oscillator (mem- 
brane clock—). At present there is a general consensus 
about the importance of Ca^"*" oscillator for SANC rate—, 
however, an important discussion still remains whether 
it is a dominant or critical factor for cardiac pacemaker 
cell functioning. Furthermore, the very existence of the 
intracellular Ca^'^ clock is not captured by the most part 
of existing essentially membrane-delimited cardiac pace- 
maker cell numerical models. 

Recently, Maltsev and Lakatta ^ have developed a new 
numerical SANC model (ML-model) featuring interac- 
tions of SR-based Ca^"*" and membrane clocks to explore 
novel mechanistic insights into cardiac impulse initia- 
tion. They started with a well-known simplified model 
of the cell structure consisting of four compartments: 
sub-sarcolemmal space (subspace), cytosol, network SR 
(nSR), and junctional, or luminal SR (jSR) (Fig. 1). As 
in most existing models the authors used an effective 
medium theory, where Ca^'^ concentrations in the sub- 
space and in jSR {C&ss and Ga.jSR) are main governing 
parameters that obey standard reaction-diffusion equa- 
tions, while RyR gating is usually considered in a simpli- 
fied manner through a dependence of the release on the 
Ca^"*" concentrations. The ML-model adopted the for- 
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FIG. 1: Fig. 1. Schematic illustration of the cell compart- 
ments and Ca^^ fluxes in the SR Ca^"*" clock toy model (see 
below Eqs.(7)-(10). SERCA is a sarco/endoplasmic reticulum 
Ca^^'^-ATPase that transfers Ca^^ from the cytosol to the SR. 



mulation of cardiac RyR function developed by Shannon 
tt al.—) and the Kurata et al. model— of primary rab- 
bit SANC. Finally the model was formulated in terms 
of a system of 29 first-order differential equations. The 
isolated SR can indeed operate as a self-sustained Ca^+ 
oscillator, described by a simple " release-pumping-delay" 
mechanism: a small spontaneous Ca^"'" release from jSR 
to the subspace occurs as the primary or initiating event. 
When Ca55 increases to a sufficient level, it amplifies 
the Ca^+ release via the mechanism of the Ca^'^^-induccd 
Ca^''" release (CICR)— ; this relatively strong, secondary 
Ca^+ release simultaneously depletes (i.e., resets) jSR. 
The released Ca'^^ is pumped into the nSR. The delay 
between releases is determined by the Ca^"'" pumping rate 
and Ca^"*" diffusion from the subspace to cytosol and also 
from nSR to jSR. As Ca^sfl slowly increases, RyRs are 
restituted, and the next release is ultimately initiated, 
etc. 

The ML-model— of coupled oscillators seems to repro- 
duce basically all recently discovered new behavioral de- 
tails of cardiac pacemaker cell function, however, this 
phenomenological integrative model ignores many impor- 
tant physiological features of the cardiac cells, in par- 
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ticular, the fine spatiotemporal structure of tfie Ca^+ 
release. Tlie model of integrated Ca^"*" dynamics does 
not describe stochastic, locally propagating Ca^"*" re- 
leases within the subsarcolemmal space. Indeed, RyRs in 
SANC, as in other cardiomyocites, seem to be arranged 
in clusters under sarcolemma (Fig. 1) and thus proba- 
bly form subsarcolemmal calcium release units (CRUs). 
In this instance, RyRs release Ca^+ into a relatively 
small volume of subspace where individual jSRs of CRUs 
approach sarcolemmai. Thus a realistic modeling of 
the Ca^"*" oscillator and SANC function should include 
a stochastic mechanism of a local Ca^"*" release gener- 
ation by CRUs. Furthermore, main assumption of the 
ML-model— , that is the Ca^"*" released from RyR chan- 
nels activates the RyR channels as like as the trans- 
sarcolemmal Ca^+ from the L-type channel in a close 
apposition, seems to be questionable. 

Recently we have applied well-known electron- 
conformational (EC) model (see, e.g.,Ref.^) for a single 
RyR channel-i2,-i^ that was shown to capture important 
features of the individual and cooperative behaviour of 
RyRs in ventricular myocytes. The EC model of RyR 
functioning under Ca^"*" stimuli is based on a biophysi- 
cal adaptation of the well-known theory of photo- induced 
structural phase transitions, which has been successfully 
applied to different solids—. Hereafter, in the Letter we 
will show that EC model of RyR channel is able to ex- 
plain and describe on equal footing a puzzling sponta- 
neous oscillatory regime of the release unit both in SANC 
under normal physiological conditions and in ventricular 
myocytes under Ca^+ SR overload. 

The ion-activated RyR channel is a giant (30 x 30 nm) 
macromolecular protein complex comprising 4 subunits 
of 565 000 Daltons each—. As other ion channels it has 
a great many of internal electron and conformational 
degrees of freedom and exhibits remarkable complexi- 
ties that need to be considered when developing realistic 
models of ion permeation. Nevertheless, until recently 
most modelling efforts for RyR channels were focused on 
a simple "hole in the wall" type model with a set of dif- 
ferent (open, closed) states. Our knowledge of molecular 
mechanisms of RyR channel functioning is limited; hence 
we are forced to start with the most general "physicists'" 
approach, which is typical for protein biophysics. Such 
an approach to the modelling of biomolecular system im- 
plies its simplifying to bare essentials with guidance from 
experimental data. 

Modelling the RyR we start with a simple and a lit- 
tle bit naive picture of the massive nanoscopic chan- 
nel like an elastic rubber tube with a varying cross- 
section governed by a conformational coordinate Q and a 
light "electronic" plug switched due to Ca^^-RyR bind- 
ing/ unbinding in the subspace—. This electronic plug in- 
teracts with the conformational coordinate and acts as a 
trigger to stimulate its change and related channel cross- 
section/conductivity. In other words, we reduce a large 
variety of RyR degrees of freedom to only two: a fast 
and a slow one, conventionally termed as electronic and 




FIG. 2: Fig. 2. Left panel: Naive "tube-and-plug" model of 
RyR channel in a closed ground state. Small letters c, o are 
used for electronically closed and open states, respectively, 
and capital letters C, O for conformationally closed and fully 
open states, respectively. Right panel: Corresponding adia- 
batic potentials in EC-model of RyR. Vertical arrows point 
to Ca^""" -induced Franck-Condon (FC) electronic transitions, 
horizontal arrow points to a non-FC tunneling transition, 
downhill arrows point to a conformational dynamics. 

conformational one, respectively. Both degrees of free- 
dom are implied to be coupled to realize an EC trans- 
formation that is the electronic control of the slow con- 
formational motion. Bearing in mind the main function 
of RyR channels, we assume only two actual electronic 
RyR states: "open" and "closed", and a single confor- 
mational degree of freedom, Q, described by a classical 
continuous variable. Figure 2 illustrates the model with 
a set of representative states of the system. 

Change in the electronic and conformational states reg- 
ulate the main RyR channel function, i.e. determines 
whether the channel is "open" and permeable for Ca^+ 
ions or "closed" and impermeable to ions. Hereafter we 
assume that the conformational variable Q specifies the 
RyR channel "cross-section" or, more precisely, a perme- 
ability for Ca^"*" , while the dichotomic electronic variable 
determines its opening and closure. This allows us to 
describe the Ca^"*" flux through the RyR as follows: 

jRyR = D{Q){CajSR-Cass), (1) 

if the channel is electronically open, and JRyR = 0, if it 
is closed. Here, the permeability coefficient D{Q) refiects 
the ease with which Ca^"*" passes through an open RyR. 
Its functional dependence on the conformational coordi- 
nate should be one of the essential model assumptions. 
D{Q) is assumed to be an increasing function of the con- 
formational coordinate, varying from zero or small leak- 
age value to some saturated value Dq (0 < D{Q) < Dq), 
when Q varies from large negative to large positive mag- 
nitudes, passing through some subconductive state at 
Q — 0. As a simplest limiting case, we may consider 
the step-like dependence D{Q < 0) = 0,D{Q > 0) = Dq. 

Hereafter, we shall assume a simple harmonic approx- 
imation for the conformational energy and use a Hookes 
harmonic law for elastic potential energy: E = ^KQ^, 
where K is the effective "elastic" constant and Q = 
relates to a base state with " unstrained tube" and a bare 
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cross-section. It is worth noting that namely the EC 
model introduces the energy to be an important factor 
of the RyR functioning. 

As a starting point of the EC model algebra we in- 
troduce a simple effective Hamiltonian for a single RyR 
channel as follows—"— 



H, = -As 



K , 



(2) 



where Sz and Sx are well-known Pauli matrices, and the 
first term describes the bare energy splitting of " up" and 
"down" (electronically "open" and "closed") states with 
an energy gap A, while the second term describes their 
mixing. It is worth noting that given h = Q we arrive 
at a classical approach with a dichotomic electronic vari- 
able. The third and fourth terms in ^ describe the 
linear and quadratic contributions to the conformational 
energy. Here, the linear term formally corresponds to the 
energy of an external conformational stress, described by 
an effective stress parameter p. The last term describes 
the EC interaction with the coupling parameter a. Here- 
after we make use of the dimensionless conformational 
variable Q\ therefore all of the model parameters (A, 
h,p,K,a) are assigned energy units. Two eigenvalues of 
our Hamiltonian 

E±{Q) ^^Q^-pQ±\ [(A - aQf + (3) 

define two branches of the adiabatic, or conformational 
potential (CP), attributed to electronically closed (i?-) 
and electronically open (-E+) states of the RyR, respec- 
tively (see Fig. 2). 

The classical dynamics of the conformational coordi- 
nate we assume to obey a conventional Langevin equa- 
tion of motion 
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dQ 



E{Q) - MTQ + R{t), 



(4) 



where first term describes a total systematic conforma- 
tional force with M being an effective RyR mass (below 
M let to be unity) , F is an effective dimensionless friction 
damping constant (relaxation rate), and R is the thermal 
fluctuation force (Gaussian- Markovian noise). The last 
two terms reflect the coupling to an external environment 
of channels. 

The thermo-activated transitions are caused by the 
thermal fluctuating force and free from the Franck- 
Condon principle. As a result of thermal transitions the 
whole system will finally relax to the thermal equilibrium. 
The temperature plays an important role in overcoming 
the energy barriers. In fact we deal with a hybrid "over" 
(thermal activation) and " through" (quantum tunneling) 
barrier transfer transitions (reactions). Quantum tun- 
neling is often addressed to specify the low-temperature 
limit of the barrier transfer transition probability that 
points to the possible way to uncover non-classical be- 
haviour for the RyR channel. We assumed the resonant 
quantum tunneling takes place between two branches of 



conformational potential starting within a tunneling zone 
5Q — e centered at the CP minimum with the probability 
obeyed the effective Gamov law 



Pt, 



Poe 



(5) 



where AQ is the width, AE is the height of the energy 
barrier, or the energy separation between the tunneling 
points and the point of the two branch intersection (see 
Fig. 2), and Pq effective tunneling attempt frequency. 

At present there is no sufficient understanding of the 
mechanisms that regulate local Ca^+ signaling in heart's 
cell despite persistent efforts to discriminate between the 
cytosolic and luminal Ca^+ activation hypotheses. Cal- 
cium enters the subspace by two pathways: across the 
sarcolemma via L-type channels and from the SR via 
RyR channels. Activation of RyRs by Ca.jSR has been 
attributed to either Ca^+ feedthrough to high-affinity cy- 
toplasmic Ca^"*" activation sites or to Ca^"*" regulatory 
sites on the luminal side of the RyR. However, most of 
experimental observations on cardiac RyRs more difficult 
to reconcile with the Ca^"*" feedthrough effect—. A lu- 
minal Ca^"*" sensor appears to continuously regulate the 
functional activity of the SR Ca^"*" stores by linking SR 
Ca^+ content to the activity of the RyRs—. However, in 
contrast with purely electronic effect of Cass, the effect 
of relatively slowly varying CajSR on the RyR channels 
is likely to be purely "mechanical" one, through the re- 
spective conformational strain applied to RyR channels. 
The effect can be naturally incorporated into EC model, 
if we assume the strain parameter p in the RyR Hamil- 
tonian to be a function of CajsR or the jSR-to-subspace 
Ca^"*" gradient. We assume p to rise with the luminal 
Ca^+ concentration in accordance with the Hill curve—: 
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-1 <p < +1 



(6) 



Ca 



where Kca is the half maximal value, n is a Hill coeffi- 
cient. Rise oi p in the interval (-1,4-1) leads to a crucial 
modification of CP from that of stabilizing closed RyR 
state to that of stabilizing open RyR state. It should 
be noted that in terms of the EC model one might in- 
troduce a critical SR load that specifies a critical ef- 
fective strain pcr = (| — -^A) when the minimum of 
the CPc branch for the electronically closed RyR state 
crosses CPq branch thus destroying the RyR bistability 
conditions and making RyRs stay close to their activation 
threshold where their cC^oO activation can presumably 
be easily provoked. 

Cooperative dynamics of the RyR clusters in CRU 
has been studied in a series of model simulations for 
11x11 square RyR lattice in our prior paper—. We 
simulated both " in vitro" dynamics of RyR lattice when 
the calcium surrounding was assumed to be simply as a 
source of effective external fields, such as effective strain, 
and "in vivo" situation when the RyR lattice dynam- 
ics is incorporated into the whole cell calcium dynamics 
(EC-CRU model). In such a case the EC dynamics of 
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RyR lattice is assumed to speciiy the number of open 
channels Nopen which in its turn specifies the release 
flux from the SR to subspace through RyR channels: 
Jrei = kreiNopeniCajSB^ - Cttss)-, where krei is & sin- 
gle RyR channel release velocity coefficient. The system 
demonstrated four different modes of behaviour, depend- 
ing on the SR load^^: inactivation mode I at a rather 
small SR load (heavily underloaded SR); single channel 
activated mode II (Ca synapse, or quark mode) at under- 
loaded SR; domino-like firing-termination mode III with 
high degree of cooperativity (cluster bomb mode) at opti- 
mally loaded SR. The cluster bomb mode is characterized 
by a step-by-step opening of the neighboring RyRs with 
formation of a cluster composed of up to 40% open RyRs. 
The process results in an effective high-gain Ca^+ release 
from the SR. However, the decrease in SR load leads to 
a lowering of the effective strain, accompanied by the 
shift of the system to preferable channel's closure. This 
negative feedback effect results, first, in a slowing down 
of the nucleation process and, second, in an evolution 
of the inverse domino-like effect with full collapse of the 
cluster of open RyRs and termination of Ca^+ release. 
Puzzlingly, the CRU simulations^ showed that SR over- 
load can result in the excitation of the RyR lattice auto- 
oscillations when the domino-like opening of a cluster of 
RyR channels resulted in effective Ca^+ release that, in 
turn, caused a lowering of the effective strain and simul- 
taneous closure of RyRs before the SR load started ris- 
ing. However, until C&jSR approached the initial value, 
close to a critical concentration all channels re-opened si- 
multaneously, and Ca^"*" release repeated spontaneously. 
This behaviour is repetitive, i.e. the system turned out 
to behave in an auto-oscillatory mode IV with a sponta- 
neous SR Ca^"*" release. It should be emphasized that in 
contrast with the ML model, the auto-oscillation mode 
IV occurs as a result of the CajSfl-dependent shape of 
the CP and a purely conformational cC -H- oO transfor- 
mation without any L-type channel activity and Ca^^- 
induced electronic coo transitions. To the best of our 
knowledge it was a first quantitative model for the CRU 
oscillatory mode. 

Obviously the abnormal Ca^^ automaticity in SR over- 
loaded ventricular myocytes suggests the normal pace- 
maker SANC activity may have some similar features and 
actually proceed with elevated Ca^+ level. Indeed, Vino- 
gradova et al.^ found a minimum diastolic Ca^^" level 
in rabbit SANC of ~200nM, that is twice the diastolic 
Ca^+ level in resting ventricular myocytes. 

Hereafter in the paper we address a more realistic EC- 
CRU model relevant for the SANCs which incorporates 
two main features of the Ca'^"'" machinery, that is a lo- 
cal and stochastic character of the Ca^+ release. As in 
our prior paper we address the CRU model with II x II 
square lattice of RyRs, however, for simplicity we neglect 
any RyR-RyR coupling. At variance witbi^ each RyR 
channel was described by a diabatic CP with non-FC tun- 
neling transition between two CP branches. The calcium 
fluxes in a simplified cell model (Fig. I) were assumed to 



obey a standard system of four differential equations—: 
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(7) 



(8) 



(9) 



(10) 



kref {CanSR — CajSR), Jr. 



where Jref 

are diffusion fluxes between the nSR and the jSR, 
between the jSR and subspace, between subspace and 
cytosol, respectively, Jup — x'^^'+ca uptake flux; 
k'p are volume ratio constants between a and /? cell 

0.12, fcP = 0.022, kl^sR = 
9.7). It should be noted that all the 



openkrei {CajSR — C'ass), Jdiff = kdiff (Cass " C'a^) 
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parameters were chosen rather as typical for integrative 
cell models than for a single CRU model. In other 
words, all the cell RyRs we consider to form a system 
of identical CRU's functioning in concert. We set initial 
conditions as follows: Ca^ = 0.1 ^M, Cass = 0.5 fiM, 
Ca,jSR = 0.05 toM, Ca„sj^ = 1.5 mM, for other Ca^+ dy- 
namics parameters:if„p = 0.0006 toM; kdiff = 5ms^^. 
As in the ML model the Cass and Cajsr time courses 
strongly depend on the refill and release rate constants, 
that are expected both to change from their original 
values of ventricular myocytes and likely mediate the 
regulation of SR Ca^^ clock ticking speed—. Main 
results of numerical simulations are presented in Fig. 3. 

As in ML model, three types of the steady-state were 
found: 1) steady rhythmic oscillations; 2) no oscilla- 
tions or damped oscillations; and 3) chaotic oscillations. 
The highest rates are reached when the oscillator ap- 
proaches to dynamic equilibrium (steady release, Fig. 3d, 
right part). The lowest rates are reached when the os- 
cillator approaches to static equilibrium. Static equilib- 
rium occurs when either the release rate is too small or 
Ca^"*" pumping rate is too fast, i.e., jSR becomes highly 
loaded. Fig. 3 distinctly shows that CRU can operate as 
a self-sustained Ca^"*" oscillator, however, the triggering 
of the Ca^"*" release is determined by the bringing Ca^sfl 
near the critical concentration, or activation threshold 
rather than by a small initial spontaneous Ca^+ release 
from jSR as in the ML-model^. Maximal Ga,jSR appears 
to slightly rise with kref, while its oscillation amplitude 
strongly depends on kref, however, in an absolutely dif- 
ferent way than in ML-model. Indeed, for a rather wide 
range of kref the Ca^si? oscillation amplitude rises with 
the refill rate, however, without full depletion of the jSR 
as it occurs in ML-model. Furthermore, in contrast with 
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FIG. 3: Fig. 3. A: Illustration of a repetitive shaping of the RyR's CP branches during refill-release mode and a flash pattern 
of open and closed RyRs in our 11x11 cluster. B: Model simulations illustrating that CRU can operate as a self-sustained 
Ca^"*" oscillator. C: SR Ca'^"^ oscillator operates following a Bowditch-like pattern: "the faster rate, the stronger release". 
Bottom panel D: Cass and Ca^sfl courses given different refill rate constant fcre/ = 10""* 5 • 10~^ ms~^ at krei ~ 0.01 ms~^ , 
kdif ~ 5ms^^, Pup — 1 fiM/ms. We make use of the EC model parameters as follows: A = 0, if = 12, a = 5 , F = 7, Ptun ~ 
0.1, Atun = 20, e = 0.001. 



the ML-model the EC-CRU model reveals the Bowditch- 
like behavior: "the faster rate, the stronger release". In- 
deed, vifithin oscillatory regime both the frequency and 
amplitude of Cass do increase vifith kref- A detailed 
analysis shows this can be explained as a result of some 
retardation effect due to a slow conformational dynam- 
ics near the CP minima; optimal condition for the inter- 
branch tunneling and Ca^"*" release are realized at higher 
kref and consequently, at higher Caj sr- Indeed, a high- 
speed increase of the Ca,jsn results in a rapid rise of 
effective stress p with a shift of the CP branches with 
its minima up and right (see Fig. 3). Given rather slow 
conformational dynamics the CP shift leads to a tunnel- 
ing retardation effect as the conformational coordinate 
Q needs a time to reach the new tunneling zone e in the 
vicinity of new CP minimum. In other words, the EC 
system is not in time for adjusting to the rapid shift of 
the CP branches. However, the higher Ca^si? the slower 
the CP shift, that allows conformational coordinate Q to 



run tunneling zone down, and provide the optimal con- 
dition for the inter-branch tunneling and Ca^+ release. 
Despite the tunneling retardation, the bigger kref the 
shorter time to reach the critical CajSR level. 

Just the opposite behavior, "the higher the rate, the 
lower the amplitude," represents the major limitation of 
the isolated Ca^+ oscillator in frames of ML-model: it is 
unable to generate high-amplitude oscillations at higher 
rates. However, the authors— have shown that the inter- 
actions of the SR and sarcolemma clocks overcomes the 
limitation of the isolated SR Ca^+ clock, i.e., the oscilla- 
tion amplitude of the full system Ca^"*" clock raises as the 
oscillation rate increases, i.e., the Bowditch phenomenon 
restores. At the same time the "Bowditch-like" behavior 
(the faster the stronger) for the both Cass and C&jSR 
appears to be a distinct feature of the EC-CRU model. In 
other words, it seems the SR Ca^+ oscillator can serve as 
a dominant source of persistent heart's beating. Thus the 
interrelation between two types of cell oscillators needs 
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in a futhcr examination. 

In summary, at variance with the integrative model 
by Maltsev and Lakatta simple biophysically based EC 
model of RyR channel describes stochastic local Ca^"*" 
releases within the subsarcolemmal space as a result of 
conformational transformations followed by a tunneling 
between two CP branches. The model is able to explain 
and describe the spontaneous oscillatory regime of the 
CRU both in pacemaker cells under normal physiological 
conditions and in ventricular myocytes under Ca^"*" SR 
overload including its subtle features such as the ampli- 
tude and frequency fluctuations. At variance with the 



integrative AlL-model, the oscillation amplitude of the 
intracellular Ca^+ clock raises as the oscillation rate in- 
creases, thus providing the Bowditch law functioning of 
the pacemaker cell without any membrane clock assis- 
tance. Despite the EC model is intentionally simplistic, 
it offers novel insight into potential mechanisms govern- 
ing by the Ca^"*" fluxes and may thus provide a starting 
point for further exploration of physical principles guid- 
ing cardiac cell functioning in vitro and in vivo. 
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